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Abstract 

In this article, we perform an systematic study of the radiative transitions among the 
bottomonium states using the heavy quarkonium effective Lagrangians, and make predictions 
for the ratios among the radiative decay widths of a special multiplet to another multiplet. 
The predictions can be confronted with the experimental data in the future. 
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1 Introduction 

In recent years, the Babar, Belle, CLEO, CDF, DO and FOCUS collaborations have discovered (or 
confirmed) a large number of charmonium-like states and revitalized the interest in the spectroscopy 
of the charmonium states [1] [2l |3l |4] . There are also some progresses in the spectroscopy of the bot- 
tomonium states. In 2004, the CLEO collaboration observed the T(1D) states in the four-photon 
decay cascade, T(3S) ^ jXb{2P), Xb{2P) ^ 7T(1D), T(1D) ^ JXb{lP), Xb{lP) ^ 7^(18), and 
obtained the mass Mx(i3D2) = (10161.1 ± 0.6 ± 1.6) MeV [S]. In 2008, the Babar collaboration 
observed the ?7&(1S) in the radiative decay T(3S) — > 7r/b(lS) fB^, and latter confirmed it in the 
radiative decay T(2S) — )■ 7?76(1S) [7]. In 2010, the Babar collaboration observed the T(l'^Dj) state 
through the decay chain T(3S) — 77T(l'^Dj) — >■ 777r+7r~T(lS) with j = 1, 2, 3, and obtained the 
mass AfT(i3D2) = (10164.5 ± 0.8 ± 0.5) MeV [8]. In 2011, the Belle collaboration reported the ob- 
servation of the spin-singlet bottomonium states /i;,(lP) and /ib(2P), which are produced in the re- 
actions e+e^ — > /if,(nP)7r+7r^ with significances of 5.5(7 and 11.2 c, respectively [5]. The measured 
masses are A//,,(ip) = (9898.25 ± 1.06ti;g?) MeV and A'4,(2P) = (10259.76 ± 0.64;j;^i^) MeV, re- 
spectively. Recently, the ATLAS collaboration observed the Xbj(3P) with j = 1, 2 in the radiative 
transitions X6j(3P) 7T(1S), 7T(2S) in the proton-proton collisions at the Large Hadron Collider 
(LHC) at the energy y/s = 7TeV HU]. The measured mass barycenter is (10530 ± 5 ± 9) MeV, 
and the hyperflne mass splitting is fixed to the theoretically predicted value of 12 MeV. And more 
bottomonium states would be observed in the future at the Tevatron, KEK-B, RHIC and LHCb. 

On the other hand, there have been several theoretical works on the spectroscopy of the bot- 
tomonium states, such as the relativized potential model (Godfrey- Isgur model) [llj . the Cor- 
nell potential model, the logarithmic potential model, the power-law potential model, the QCD- 
motivated potential model |12j, the relativistic quark model based on a quasipotential approach 
in QCD |13| . the Cornell potential model combined with heavy quark mass expansion ^14) . the 
screened potential model [15], the potential non-relativistic QCD model [16], the confining poten- 
tial model with the Bethe-Salpeter equation [T7|, etc. In Table 1, we list the experimental values 
of the bottomonium states compared with some theoretical predictions [9l [TOl [TT] [151 HI] • 

The charmonium and bottomonium states have analogous properties, the hadronic transitions 
and radiative transitions among the heavy quarkonium states have been studied by the QCD 
multipole expansion [T^ , the nonrelativistic potential model [501 HIl 1221 [23] , the heavy quarkonium 
effective theory [24], the coupled-channel approach [ST] [22], the hybrid approach based on the 
multipole expansion and heavy quark symmetry |25j . etc. In the nonrelativistic potential models, 
the El and Mi transitions among the bottomonium states are usually studied by the following 
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where the E-y is the photon energy, the Ef is the energy of the final bottomonium state, the Mi is 
the mass of the initial bottomonium state, and the angular matrix factor Cfi is 



The values of the matrix elements (n'^'* ^^^'j' I I n^^^^Lj) and {n'"^" ^^L'j' \ ii?^^^Lj) and their 
/<? corrections depend on the details of the wave-functions, which are evaluated using a special 
potential model [Til [HI HSl [E] ■ In Ref. [26] , we focus on the traditional charmonium scenario of the 
new charmonium-like states and study the radiative transitions among the charmonium states with 
the heavy quarkonium (or meson) effective theory based on the heavy quark symmetry [24[ [?7l [5S] , 
which have been applied to identify the excited Dg and D mesons, such as the 1)5(3040), -Ds(2700), 
1)5(2860), I?(2550), ^(2600), ^(2750) and £'(2760) [HI 13011311131 133] ■ In this article, we extend 
our previous works to study the radiative transitions among the bottomonium states using the 
heavy quarkonium effective theory. 

The article is arranged as follows: we study the radiative transitions among the bottomo- 
nium states with the heavy quarkonium effective Lagrangians in Sect. 2; in Sect. 3, we present the 
numerical results and discussions; and Sect. 4 is reserved for our conclusions. 

2 The radiative transitions with the heavy quarkonium ef- 
fective Lagrangians 

In the infinite heavy quark mass limit, the heavy quarkonium states do not have heavy quark fiavor 
symmetry and spin symmetry, while for the intermediate heavy quark mass, the heavy quark spin 
symmetry is expected to make sense [34l[35]. In fact, the c and h quarks have large but finite masses, 
we can construct heavy quarkonium effective Lagrangians based on the heavy quark spin symmetry. 
In the nonrelativistic QCD for the heavy quark systems, we introduce three typical energy scales 
™Qj "^Q^^j where ttiq and v are the heavy quark masses and velocities respectively, and 

count the operators with the power orders of w, take the heavy quarkonia as bound states and study 
them with the nonrelativistic Schrodinger field theory, and apparent Lorentz covariance is lost. In 
the heavy quark effective theories, we introduce two typical energy scales mq and Aqcd, count the 
operators with the power orders of l/mg, and take the heavy quarkonia as the basic relativistic 
quantum fields rather than bound states at the hadron level, the calculations are more simple, and 
apparent Lorentz covariance is maintained. In the two approaches, the heavy quarkonium states 
are classified in the same scheme. And the two approaches both have advantages and shortcomings. 

The bottomonium states can be classified according to the notation n^*+^Lj, where the n is 
the radial quantum number, the L is the orbital angular momentum, the s is the spin, and the 
j is the total angular momentum. They have the parity and charge conjugation P — (— 1)^+-'^ 
and C = (—1)^+'*, respectively. The states have the same radial quantum number n and orbital 
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States 


Experimental data jgi[TUl[T5] 


Theoretical values [TBI 


Theoretical values [llj 


IS 


T(l^Si) 


9460.30 ± 0.26 


9460 


9460 




'7b(l^So) 


9390.9 ± 2.8 


9389 


9400 


2S 


T(2^Si) 


10023.26 ± 0.31 


10016 


10000 




^^(2'So) 




9987 


9980 


3S 


T(3^Si) 
r?b(3iSo) 


10355.2 ±0.5 


10351 


10350 






10330 


10340 


4S 


T(4^Si) 


10579.4 ± 1.2 


10611 


10630 






10595 




5S 


T(5^Si) 


10865 ± 8 


10831 


10880 






10817 




6S 


T(6^Si) 
'7b(6^So) 


11019 ± 8 


11023 


11100 






11011 




7S 


T(7^Si) 




11193 








11183 




IP 


X62(1^P2) 


9912.21 ±0.26 ±0.31 


9918 


9900 




X&i(l'Pi) 


9892.78 ±0.26 ±0.31 


9897 


9880 




Xw(l'Po) 
/ifc(liPi) 


9859.44 ± 0.42 ± 0.31 


9865 


9850 




9898.25 ± 1.06l};^? 


9903 


9880 


2P 


X62(2^P2) 


10268.65 ±0.22 ±0.50 


10269 


10260 




X6i(23Pi) 


10255.46 ±0.22 ±0.50 


10251 


10250 




Xbo(23Po) 
/ib(2iPi) 


10232.5 ± 0.4 ± 0.5 


10226 


10230 




10259.76 ±0.64li;oi^ 


10256 


10250 


3P 


X62(3^P2) 


10536±5±9 


10540 






Xw(33Pi) 


10524±5±9 


10524 






Xbo(33Po) 




10502 
10529 




4P 


Xb2(4^P2) 

X6i(43Pi) 
X6o(43Po) 
/ifc(4iPi) 




10767 
10753 
10732 
10757 




5P 


X62(5^P2) 

X6i(53Pi) 
Xbo(53Po) 
/ih(5iPi) 




10965 
10951 
10933 
10955 




ID 


T3(1^D3) 




10156 


10160 




Tzll^Da) 


10161 ± 0.6 ± 1.6 


10151 


10150 




T(13Di) 




10145 


10140 




7762(1^D2) 




10152 


10150 


2D 


T3(2^D3) 




10442 


10450 




T2(23D2) 




10438 


10450 




T(23Di) 




10432 


10440 




?762(2lD2) 




10439 


10450 


3D 


T3(3^D3) 
T2(3'^D2) 




10680 
10676 






T(33Di) 

??62(3lD2) 




lUD / U 








10677 




4D 


T3(4^D3) 

T2(43D2) 

T(43Di) 

??62(4lD2) 




10886 
10882 
10877 








10883 




5D 


T3(5^D3) 

T2(53D2) 




11069 
11065 






T(53Di) 
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momentum L can be expressed by the superfields J, J^, J'^", etc 
1 + ii^ 1 - y 
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where the denotes the four velocity associated to the superfields. We multiply the bottomonium 



fields T3^„, 



Y2 y 2 



with the factors ^yMrl 
3 



/Mr, ■ ■ • ' 

respectively, and they have dimension of mass |. The superfields J, J^, J^'' are functions of the 
radial quantum numbers n, the fields in a definite superfield have the same n, and form a multiplet. 
The superfields JA'i - a'i. have the following properties under the parity, charge conjunction, heavy 
quark spin transformations. 
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where S*, S" € SU{2) heavy quark spin symmetry groups, and [S, f] = [5", /] — 0. The S'-wave multi- 
plet contains the spin-singlet bottomonium states; while the P-wave and D-wave multiplets contain 
both the spin-singlet and spin-triplet bottomonium states, there exist mass splittings among (or 
between) the bottomonium states in the same multiplet. In the nonrelativistic quark models, we 
resort to the spin-spin, spin-orbit and tensor interactions to take into account the mass splittings. 
In the present case, we can use the superfields J, J'^, J^'" and the Dirac matrix a^i, to construct 
the heavy quarkonium effective Lagrangians at the hadronic level to account for the mass-splittings 
in a multiplet [35] . The heavy quark effective Lagrangian can be written as 



hyiv • Dhy 



1 



2mQ 



(5) 



where = — v^'-v ■ D, the Z?^ is the covariant derivative, and the hy is the heavy quark 
field. The heavy quark spin symmetry breaking terms are of the order 0{l/mQ), Oil/mQ), • • • , 
and the corresponding corrections to the heavy quarkonium states (or mesons) can be counted as 
^QCo/'mQ, {^QC D / fnof' , • • • , respectively, here we introduce the scale Aqcd to characterize the 
bound states. The mass-splittings in a multiplet can also be reproduced with the heavy quarkonium 
effective Lagrangians at the next-to-leading order [35]. 

The radiative transitions between the m and n bottomonium states can be described by the 



4 



following heavy quarkonium effective Lagrangians [Ml UHl 136^1 , 

Css = ^^(m,n)Tr[J(m)a^,J(n)]F'^^ 

m.'fi 

CsP = ")Tr [J(m) J^(n) + J^(n) J(m)] , 

CPD = ^5(m,n)Tr[V(m)J^(n) + 7=^(71) J^,(m)] l^'', (6) 

where Jfj,^...f^^ — J^Jl^...i_i^l° , V'^ — F^^Vi^, the F^" is the electromagnetic field tensor, and the 
S{m,n) are the coupling constants and have different values in the £ss, ^SP and £pD- 

The Lagrangians Csp and CpD preserve parity, charge conjugation, gauge invariance and heavy 
quark spin symmetry, while the Lagrangian Css violates the heavy quark spin symmetry. The 
effective Lagrangians Csp and Cpd describing the electric dipole Ei transitions can be realized in 
the leading order 0(1) in the heavy quark effective theory, while the effective Lagrangian Css de- 
scribing the magnetic dipole Mi transitions can be realized in the next-to-leading order O^I/ttiq). 
The corrections to the Csp, Cpjj and Css come from the next-to-leading order ©(l/mg) and 
the next-to- next-to-leading order ©(I/ttiq), respectively. We can construct the corresponding La- 
grangians with the superfields J, J^, J^'^ , the Dirac matrix ct^^, the electromagnetic field tensor 
Fpjy and the four- vector v^, and introduce additional unknown coupling constants, which have to 
be fitted to the precise experimental data in the future, and study the spin symmetry violations 
in the radiative decays to the bottomonium states in an special multiplet. 

In the case of the charmonium states, the 0{v^ /(?) corrections to the Ei transitions come from 
the magnetic quadrupole M2 and electric octupole E^ terms, the current average values of the 
relative amplitudes are M^/^jEf + M| + El = (-10.0 ± 1.5) x 10"^ and E3/ y/ Ef + + E^ = 
(1.6 ± 1.3) X 10-2 in the Xc2(lP) ^ 7J/i^ decay, M2 / ^/W^^Mf+E| ^ (1.0 ± 1.4) x lO^^ and 
Es,/^Ei + M'i + El ^ (-1.0±l.l)xlO-2 in the i/;(2S) ^ 7Xc2(lP) decay, M2/ E'f + M'i + El = 
{-bAt\l) X 10-2 in the Xci(lP) ^ iJ/i^ decay, M2 / V'^Ef+MfTsf = (2.9 ± 0.8) x 10^^ in 
the ^"(28) — > 7Xci(lP) decay from the Particle Data Group [T^. Although the values differ from 
the theoretical expectations w^/c^ « 0.3, the corrections are very small, we expect that the corre- 
sponding corrections for the bottomonium states are also very small. 

In the screened potential model 15 , the predictions for the decay widths of the transitions 
T(2S) — > Xfcj(lP)7i T(3S) Xfcj(2P)7, j — 0,1,2, become better compared to the experimental 
data after the 0{v'^ /c^) corrections are taken into account, the 0{v'^ /c^) corrections are not large; 
on the other hand, the 0{v'^ /c^) corrections in the radiative transitions 4S — 2P, 4S — IP, • • • , 
are very large, and it is not a good approximation for taking the spin symmetry breaking effects 
perturbatively. At the hadron level, large spin symmetry violations mean that the ©(I/toq) 
corrections are large enough to ruin the leading order approximation, while we do not have enough 
experimental data to fit the unknown coupling constants in the phenomenological Lagrangians, and 
those parameters cannot be canceled out with each other to result in independent ratios among the 
radiative decay widths; we expect that neglecting the next-to-leading order and next-to-next-to- 
leading order corrections for the Csp, CpD and Css respectively leads to uncertainties of the order 
0{Kqcd /fTib)] in other words, we expect that the flavor and spin symmetry breaking corrections 
of the order ©(l/mg) to the effective Lagrangians Csp and Cpo are smaller than (or not as large 
as) the leading order contributions. 

In the heavy quark limit, the contributions of the order ©(l/mg) are greatly suppressed. 
For example, the branching ratios of the radiative decay widths of the T(2S) have the hierar- 
chy Br(T(2S) Xfcj(lP)7) 10"^ > Br (T(2S) -> 776(18)7) - 10-^ j = 0,1,2 [18], although 

^The Lagrangian £p£) was introduced by F. De Fazio in Ref. 1361 for the first time, I failed to take notice 
of the article, and constructed the Lagrangian Cpn independently in Ref. I26|. furthermore, I constructed spin 
violation Lagrangian Css ™ Ref. |26|. When the present article was submitted to |http: / / arxiv.org| and appeared as 
larXiv:1101.2506. Dr. F. De Fazio drew my attention to Ref. [36]. 
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the radiative decays T(2S) Xbj{^P)l are suppressed in the phase-space. In the case of the 
charmonium states, the decay widths r(-0' — Xcol) = 29.2 KeV, r(-0' — Xcij) = 28.0 KeV, 
r(V'' ^ Xc27) = 26.6 KeV, r(xco ^ J/V'7) = 119.5 KeV, r(xci ^ = 295.8 KeV, T{xc2 ^ 

J/ipl) = 384.2 KeV in the E'l transitions are very large compared with the decay widths — 
777c) = 1.03 KeV and r{J/ip 77/c) = 1.58 KeV in the Ah transitions [18], as the Mi tran- 
sitions take place at the order ©(I/ttiq). However, it is difhcult to count for that the decay 
widths r('!/'' — >■ 7??c) = 1.03 KeV and Tl J/ip jr]c) = 1.58 KeV are of the same order by 
the 0{l/mQ) depression. In the nonrelativistic potential models, the overlap matrix elements 

(n'^* ^^S'j' I n'^'^^^Sj) in the leading order approximation are 1 and for the transitions J/tp 777c 
and ip' ~^ I'Hc respectively, the perturbative as corrections and relativistic jc? corrections cannot 
count for the decay widths V{'^' — )• 77^0) = 1.03 KeV and T[J/ijj — >• 777c) = 1.58 KeV simultane- 
ously [IH] , we have to suppose large high order corrections or introduce new mechanisms [31 [37] . 
In the heavy quarkonium effective theory, the heavy quark spin symmetry cannot count for the 
analogous decay widths as the ip' and J /il) have the same quantum numbers except for the radial 
numbers and masses, we can use the coupling constants (5(2, 1) and (5(1, 1) to parameterize all those 
corrections in the nonrelativistic potential models, and take them as free parameters fitted to the 
experimental data, then use those parameters to study other processes. 

From the heavy quarkonium effective Lagrangians Css^ ^sp and Cpo, we can obtain the 
radiative decay widths F, 

2j + 1^87rM2l I ' 

where the T denotes the scattering amplitude, the Pcm (or fc-y) is the momentum of the final states 
in the center of mass coordinate, the J2 denotes the sum of all the polarization vectors, the j is 
the total angular momentum of the initial state, and the M is the mass of the initial state. For 
example, in the radiative decays T3(77i'^D3) Xf)2("-'^P2)7, 

where the e^^Q(A, q), e^y(A,p) and e^(A, k) are the polarization vectors of the bottomonium states 
T3(77i^D3), x&2('T-^P2) and the photon, respectively. The summation of the polarization vectors of 
the total angular momentum j = 1,2,3 states results in the following three formulae, 

A ^ 

E* _ 9i^agvi3 + g^pgi^a _ gp,uga0 

^^.u^aP — 2 3 ' 

A 

f-liup^a0T = g {gpagupgpT + gpagvTgpp + g^pguagpT + g^pgiyrgpa + g^Tguagp^ + gprgiy^gpa) 

1,^^^ ^ ^ ^ ^ ^ 

-j^ ygpagupgiiT + g^pgi^pgar + gprgupgap + guag^pgisr + guisg^pgar + gvTgppgais 

+gpagpugpr + gp/sgtu^gar + gprgpiyga/s) , (9) 

we use the FeynCalc to carry out the contractions of the Lorentz indexes. 

3 Numerical Results 

In calculations, the masses of the bottomonium states are taken as the experimental values from 
the Belle collaboration [9], the ATLAS collaboration \T0\ and the Particle Data Group [18], see 
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Table 1; for the unobserved bottomonium states, we take the values from the screened potential 
model as the physical masses [15] . For the bottomonium states above the BB threshold, the masses 
receive contributions from the intermediate meson-loops, such as the BB, BB*, B*B, B*B* , ■ ■ ■ . 
If the coupled channel effects are large enough to distort the bb configurations and induce some 
meson-meson components in the wave-functions, the mass-shifts are very large. For example, in the 
case of the charmonium states, the mass-shifts originate from the coupled channel effects are about 
— (400 — 500) MeV ^38 , we have to take the masses from the special potential quark model as the 
bare masses, and redefine the bare masses to reproduce the physical masses or the experimental 
values, the net coupled channel effects lead to the mass-shifts of the order 10 MeV. For the 
observed bottomonium states, the masses from the screened potential model [15j are consistent 
with the experimental data from the Belle collaboration [3], the ATLAS collaboration [TU] and 
the Particle Data Group [TF. Although the net coupled channel effects result in mass-shifts for 
all the bottomonium states, we expect the mass-shifts are not large enough to destroy counting 
them as small quantities, and only modify the radiative decay widths through the phase-factor A;^ 
perturbatively. Neglecting the coupled channel effects can lead to the uncertainties about Afc^/A;^. 

The numerical values of the radiative decay widths are presented in Tables 2-6, where we retain 
the unknown coupling constants 5{m,n) among the multiplets of the radial quantum numbers 
m and n. In general, we expect to fit the parameters 5{m, n) to the precise experimental data, 
however, the experimental data are far from enough in the present time. In Tables 2, 7-10, we 
present the ratios among the radiative transitions among the bottomonium states. 

The CLEO collaboration have observed that the doubly radiative decay ip' — 77 J/-0 takes place 
through the decay cascade ■0' — > JXcj, Xcj ~^ i = 0, 1, 2 with additional tiny non- resonance 

contributions [39l |40] . Recently, the BESIII collaboration observed the first evidence for the direct 
two-photon transition ip' — )■ "fjJ/ip with the branching fraction (3.3 ± 0.6^^^) x 10""* in a sample 
of 106 miUion V' decays collected by the BESIII detector [H]. In Ref . [42] . He et al study the 
discrete contributions to decay ip' — 77 J/'0 due to the Ei transitions using the heavy quarkonium 
effective Lagrangian [53] . We expect that the corresponding doubly radiative decay T(2S) — ?> 
77T(1S) occurs through the analogous decay cascade T(2S) — >■ 7x6^ (IP), Xfci(lP) ^ 7T(1S). 
Experimentally, the doubly radiative decays T(3S) — > 77T(2S) and T2(1D) — > 77T(1S) have been 
observed |18j . Once the coupling constants S{m, n) in the heavy quarkonium effective Lagrangians 
^SSj ^SP smd CpD are fitted to the precise experimental data, we can use them to study the doubly 
radiative decays, T(3S) — 77T(2S), T2(1D) 77T(1S), and other physical processes have the 
radiative transitions as their sub-processes, or study the singly radiative decays. For example, the 
radiative decays T(2S) — 7Xf,j(lP) and T(3S) — )■ 7Xbj(2P) have been observed experimentally, we 
can use them to fit the coupling constants S{2, 1) and S{3, 2), and make predictions for the decay 
widths r(77fc(2S) ^ 7/i6(lP)) and r(r;b(3S) ^ 7/ifc(2P)), see Tables 3,7. 

The widths of the radiative transitions of the 5- wave to the P-wave bottomonium states listed in 
the Review of Particle Physics are presented in Table 3 :18_. From those radiative decay widths, we 
can obtain the ratios among the radiative decay widths of the iS-wave to the P-wave bottomonium 
states, which are presented in Table 7. From the table, we can see that the agreements between 
the experimental data and the theoretical predictions are rather good, and the heavy quarkonium 
effective theory in the leading order approximation works rather well. The ratios presented in 
Tables 2, 7-10 can be confronted with the experimental data in the future at the Tevatron, KEK- 
B, RHIC and LHCb. In Tables 7-10, we also present values come from the screened potential 
model for comparison [15] . and no definite conclusion can be made. 

In calculations, we observe that the radiative decay widths are sensitive to the masses of the 
bottomonium states in some channels. The experimental value of the T(4S) from the Particle 
Data Group is il/x(4S) = (10579.4 ± 1.2) MeV ,18], while the predictions of the 4S states from the 
screened potential model are Mx(4s) = 10661 MeV and M^^(4g-) = 10595 MeV [TS]. In this article, 
we take the values Mx(4s) = (10579.4±1.2) MeV and M^^(4g) = 10595 MeV, i.e. Mx(4s) < M^^i.(4s), 
the prediction of the ratio r{Y'{4^)^xb2^{3pyl) ~ 6.297 seems rather exotic. Naively, we expect that 
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the spin-1 states T(nS) have larger masses than the corresponding ones of the spin-0 states r/b(nS). 



If we take Mx(4S) > M, 



then the ratio r(r),(4S)^/t,(3P)7) 

tnen tne ratio r(T(4S)^x62(3P)7) 



3 seems rather natural, more 



experimental data are still needed to make better predictions. 

The radiative decay widths F oc A;^ , the uncertainties originate from the masses of the bottomo- 
nium states can be estimated as 



Ar 



< 



Afc3 



M- 
IT: 




( ^M, 




/AM 



(10) 



where we have neglected the terms 



AMj 
Mi 



and 



AMf 



not enhanced by the factor 



the 



subscripts i and / denote the initial and final bottomonium states, respectively. The uncertainties 
of the masses of the bottomonium states T(5'^Si) and T(6'^Si) are 8MeV from the Particle Data 
Group [18], and the corresponding relative uncertainties are 0.074% and 0.073%, respectively, 
and expected to result in the largest uncertainties for the decay widths. The ATLAS collaboration 
fixed the hyperfine mass splitting between the Xfc2(3P) and Xbi(3P) to the theoretically predicted 
value of 12 MeV, which maybe result in larger uncertainty for the mass barycenter, as the theoretical 
and experimental values of the hyperfine mass splitting always have differences [lOj . On the other 
hand, the uncertainties of the masses of other bottomonium states are very small, about (or less 
than) 1 MeV, from the recent Belle data [9 and the Review of Particle Physics |l8j, the relative 
uncertainties are tiny. If the factor ^j^^'i is not large enough (i.e. the difference between 

the radial quantum numbers of the initial (m) and final (n) states is larger than 1, |to — n| > 1), the 
uncertainties originate from the masses of the bottomonium states are of a few percents, and can 
be neglected. For example, the uncertainties in the radiative decays of the S'-wave to the S-wave 
and the S'-wave to the P-wave bottomonium states are 

AF 
AF 
AF 
AF 
AF 
AF 
AF 



< 0.9% 


for 


T(43Si)^776(l'So)7, 


< 0.6% 


for 


T(43Si)^776(2'So)7, 


< 1.5% 


for 


T(43Si)-^%(3iSo)7, 


< 1.9% 


for 


T(53Si)^r,b(liSo)7, 


< 2.9% 


for 


T(53Si)-^r,b(2iSo)7, 


< 4.6% 


for 


T(53Si)^77b(3iSo)7, 


< 9.0% 


for 


T(53Si)->776(4'So)7, 



(11) 
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and 




— < 0.6% 



for 



T(43Si)^Xw(l'Pi)7, 



for 



T(43Si)^Xw(2'Pi)7, 




— < 6.5% 



for 



T(43Si)->Xbi(3'Pi)7, 



for 



T(53Si)^Xbi(l'Pi)7, 



for 



T(53Si)^Xw(2'Pi)7, 



for 



T(53Si)^Xbi(33Pi)7, 



for 



T(53Si)^Xw(4'Pi)7, 



(12) 



respectively. From above equations, we can see that the relative uncertainties of the decay widths 
of the T(5'^Si) — r/b(n^So), Xbi('^^Pi) ^"^^ larger than the corresponding ones of the T(4^Si) — )■ 
r]t{n^So), Xf)i(^'^Pi): as the mass of the T(5'^Si) has larger uncertainty [H]. In calculations, we 
observe that the upper bound of the uncertainties originate from the masses of the T(4'^Si) and 
T(5^Si) are about ten (or twenty) percent for the values |to — rt| < 1. 

In this article, we have neglected the corrections of the order ©(l/mg) and 0(l/mQ) for the 
heavy quarkonium effective Lagrangians Css (as the spin-flipped Lagrangian Css is of the order 
©(I/toq)) and Csp, Cpo, respectively, which maybe result in rather large uncertainties. If those 
corrections can be counted as ^.qcd I'mh-, taking the Lagrangians £sS: ^sp and Lpo to fit the 
experimental data to determine the unknown coupling constants can result in uncertainties of the 
order 0{Aqcd/ fnt). Such crude estimation maybe not work well, we should bear in mind that 
the magnitudes of the contributions from the higher order terms in the heavy quarkonium effective 
theory be determined experimentally. 

4 Conclusion 

In this article, we extend our previous work on the radiative transitions among the charmonium 
states to study the radiative transitions among the bottomonium states in an systematic way based 
on the heavy quarkonium effective theory, and make predictions for ratios among the radiative 
decay widths of a special multiplet to another multiplet, where the unknown couple constants 
6{m, n) are canceled out with each other. The predictions can be confronted with the experimental 
data in the future at the Tevatron, KEK-B, RHIC and LHCb. 
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r(T — > vh^] 




r(T^»jt,7) 


2S^ IS 


0.091 


0.163 


1.781 


3S^ IS 


0.299 


0.674 


2.254 


38 -^2S 


0.019 


0.034 


1.761 


48^ IS 


0.532 


1.408 


2.648 


4S^ 2S 


0.076 


0.207 


2.711 


4S ^ 3S 


0. ()()() 


0.017 


2.673 


5S ^ IS 


0.952 


2.290 


2.406 


58^ 2S 


0.233 


0.527 


2.261 


58-)- 38 


0.057 


0.113 


1.962 


58^48 


0.008 


0.016 


2.059 


68-> IS 


1.240 


3.277 


2.643 


68^28 


0.366 


0.973 


2.658 


68^38 


0.118 


0.308 


2.607 


68-^48 


0.029 


0.093 


3.158 


68-^ 58 


0.003 


0.004 


1.147 


78^ IS 


1.620 


4.330 


2.673 


7S^ 2S 


0.563 


1.517 


2.697 


7S^ 38 


0.224 


0.597 


2.669 


78^48 


0.079 


0.244 


3.081 


78^ 58 


0.021 


0.038 


1.839 


78^68 


0.002 


0.005 


2.204 



Table 2: The ratios of the radiative decay widths of the 5- wave to the 5- wave bottomonium 
states, where the unit of the widths is 5'^{m,n). 
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r 


X — >■ Yo'y 


T — ^ Yi T 


T — >■ Vn^ 

AU / 


Tih — >■ hh^ 


2S IP 


2.355 


2.281 


1.492 


2.176 


PC sip 


9 087 _i_ n 91 Q 


9 9n7 -1- n 999 


1 91 c; -1- n 1 R9 




QQ sip 


lOo.U 


yo. ( u 


oo. ID 




3S — > IP 


< 0.386 


< 0.035 


0.061 




3S — >^ 2P 


1.123 


1.028 


0.634 


1.084 


3S-S^ 2P 


2.662 ± 0.406 


2.560 ±0.337 


1.199 ±0.164 




4S^ IP 


447.6 


291.6 


110.9 


909.9 


4S 2P 


49.28 


33.40 


13.59 


110.7 


4S 3P 


0.143 


0.178 


0.161 


0.901 


5S^ IP 


1223 


777.3 


283.1 


1983 


5S-^2P 


326.4 


208.6 


77.14 


483.1 


5S-^3P 


58.35 


38.88 


15.55 


71.10 


5S ^ 4P 


1.627 


1.453 


0.807 


0.678 


6S IP 


1854 


1169 


419.5 


3375 


6S ^ 2P 


628.4 


396.5 


143.3 


1133 


6S-^3P 


178.4 


114.9 


43.38 


318.9 


6S^4P 


26.72 


18.80 


7.829 


49.22 


6S ^ 5P 


0.275 


0.329 


0.221 


0.552 


7S ^ IP 


2770 


173G 


()14.r) 


5003 


7S ^ 2P 


1132 


707.3 


251.2 


2025 


7S 3P 


432.2 


273.2 


99.66 


766.9 


7S-^4P 


124.2 


81.88 


31.21 


223.4 


7S-^5P 


19.91 


14.24 


5.863 


35.84 



Tabic 3: The radiative decay widths of the S-wavc to the P-wavc bottomoniiim states, where the 
unit is lQ~^5^(rn,n). The wide-hat denotes the experimental values, where the unit is KeV. 



r 


X2 T7 


Xi T7 


Xo T7 


hb Vbl 


IP^ IS 


87.27 


76.89 


60.89 


121.6 


2P^ IS 


458.6 


438.5 


402.5 


559.5 


2P-s^ 2S 


14.76 


12.55 


9.232 


20.13 


3P-> IS 


1017 


988.1 


927.5 


1181 


3P->2S 


126.5 


118.2 


103.7 


148.2 


3P-5> 3S 


6.006 


4.902 


3.241 


7.973 


4P ^ IS 


1733 


1688 


1602 


1940 


4P ^ 2S 


366.3 


347.4 


319.0 


403.1 


4P^3S 


67.28 


60.86 


51.94 


74.70 


4P^4S 


6.706 


5.331 


3.638 


4.343 


5P^ IS 


2536 


2484 


2379 


2787 


5P ^ 2S 


711.8 


683.7 


644. 7 


766.1 


5P ^ 3S 


209.0 


195.7 


178.9 


224.0 


5P -> 4S 


55.67 


50.00 


43.23 


45.56 


5P^ 5S 


1.037 


0.662 


0.328 


2.702 



Table 4: The radiative decay widths of the P-wave to the S'-wave bottomonium states, where the 
unit is 10~^(5^(m, n). 
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r 


X2 Tf 37 


X2 T27 


X2 T7 


Xi ^^27 


Xi '^1 


Xo T7 


hb 1121 


2P ID 


2.066 


0.321 


0.032 


1.089 


0.580 


1.160 


2.155 


3P ID 


74.46 


12.79 


0.963 


58.12 


21.98 


73.83 


86.58 


3P ^ 2D 


1.207 


0.244 


0.019 


0.827 


0.337 


0.597 


1.2G2 


4P ID 


293.8 


51.27 


3.676 


239.4 


86.16 


311.1 


339.6 


4P 2D 


47.27 


8.752 


0.615 


38.51 


13.57 


44.49 


52.79 


4P ^ 3D 


0.959 


0.196 


0.016 


0.595 


0.248 


0.415 


0.889 


5P ^ ID 


653.1 


114.8 


8.061 


544.6 


192.2 


719.2 


759.9 


5P -> 2D 


188.5 


34.43 


2.368 


159.1 


54.85 


197.9 


215.7 


5P -s- 3D 


32.22 


5.995 


0.424 


25.90 


9.198 


30.28 


35.64 


5P -s- 4D 


0.719 


0.149 


0.012 


0.429 


0.176 


0.307 


0.649 



Tabic 5: The radiative decay widths of the P-wave to the D-wave bottomonium states, where the 
unit is 10~'*(5^(m, n). 



r 


X27 


T2 X27 


T2 Xi7 


^ X27 


Xi7 


T Xo7 


112 hi 


ID -)> IP 


14.48 


3.847 


14.38 


0.351 


6.664 


12.78 


16.28 


2D IP 


138.9 


34.00 


113.0 


3.638 


60.85 


96.19 


147.2 


2D ^ 2P 


5.302 


1.237 


4.632 


0.123 


2.331 


4.460 


5.852 


3D ^ IP 


400.7 


98.87 


317.3 


10.64 


172.7 


258.3 


416.3 


3D 2P 


67.05 


16.30 


53.59 


1.730 


28.56 


44.49 


69.83 


3D ^ 3P 


3.063 


0.704 


2.697 


0.069 


1.330 


2.687 


3.322 


4D IP 


781.8 


193.7 


610.7 


20.89 


335.2 


487.6 


804.2 


4D ^ 2P 


216.5 


53.17 


169.2 


5.737 


91.91 


135.6 


222.2 


4D ^3P 


41.97 


10.15 


33.61 


1.079 


17.92 


28.49 


43.37 


4D 4P 


1.740 


0.393 


1.658 


0.038 


0.819 


1.738 


2.062 


5D IP 


1261 


313.5 


977.0 


33.69 


537.8 


769.1 


1289 


5D ^ 2P 


453.4 


111.9 


350.5 


12.10 


191.6 


275.G 


462.0 


5D ^ 3P 


142.4 


34.85 


111.4 


3.751 


60.26 


90.07 


145.3 


5D ->■ 4P 


27.29 


6.563 


22.51 


0.693 


11.93 


19.30 


29.18 


5D -S- 5P 


1.166 


0.259 


1.149 


0.025 


0.559 


1.173 


1.415 



Table 6: The radiative decay widths of the D-wa.ve to the P-wave bottomonium states, where the 
unit is 10~^^^(m, n). 
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r 








2S ^ IP 


0.969 [0.846] 


0.634 [0.451] 


0.924 [2.264] 


2S ^ IP 


0.965 ±0.134 


0.531 ±0.087 




3S-^ IP 


0.679 [0.111] 


0.277 [0.040] 


1.670 [4.508] 


as ^ 2P 


0.915 [0.803] 


0.565 [0.405] 


0.966 [3.322] 


3S ^ 2P 


0.962 ±0.194 


0.450 ±0.092 




4S ^ IP 


0.651 [0.233] 


0.248 [0.001] 


2.033 [6.558] 


4S ^ 2P 


0.678 [0.002] 


0.276 [0.375] 


2.246 [3.857] 


4S ^ 3P 


1.245 [1.423] 


1.125 [1.038] 


6.297 [24.81] 


5S ^ IP 


0.636 


0.231 


1.622 


5S ^ 2P 


0.639 


0.236 


1.480 


5S ^ 3P 


0.666 


0.266 


1.218 


5S ^ 4P 


0.893 


0.496 


0.417 


6S ^ IP 


0.631 


0.226 


1.820 


6S ^ 2P 


0.631 


0.228 


1.803 


6S ^ 3P 


0.644 


0.243 


1.788 


6S ^ 4P 


0.703 


0.293 


1.842 


6S ^ 5P 


1.194 


0.802 


2.007 


7S IP 


0.627 


0.222 


1.806 


7S ^ 2P 


0.625 


0.222 


1.790 


7S ^ 3P 


0.632 


0.231 


1.775 


7S ^ 4P 


0.659 


0.251 


1.799 


7S ^ 5P 


0.715 


0.294 


1.800 



Table 7: The ratios among the radiative decay widths of the S'-wave to the P-wave bottomo- 

nium states, where F = ^ r, F(T — > Y27) — ^'^Z-^'^^'^l — 1, and the wide-hat denotes the 

experimental values, the values in the bracket come from the screened potential model |15| . 
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1 1 


1 1 


"■6 ^ '/67 




n SSI fn 0901 




1 ^Q4 f1 11 41 


OP _v 1 S 




n S7S fn "^finl 


1 990 fl 44nl 


2P ^ 2S 


0.850 [0.972] 


0.625 [0.817] 


1.364 [1.077] 

L J 


3P -> IS 


0.972 [0.501] 


0.912 [0.127] 


1.161 [1.399] 


3P ^ 2S 


0.934 [0.782] 


0.820 [0.533] 


1.172 [1.495] 


3P ^ 3S 


0.816 [0.898] 


0.540 [0.691] 


1.327 [1.045] 


4P IS 


0.974 


0.924 


1.120 


4P^ 2S 


0.948 


0.871 


1.101 


4P 3S 


0.905 


0.772 


1.110 


4P ^ 4S 


0.795 


0.543 


0.648 


5P ^ IS 


0.979 


0.938 


1.099 


5P -> 2S 


0.960 


0.906 


1.076 


5P ^ 3S 


0.936 


0.856 


1.072 


5P ^ 4S 


0.898 


0.776 


0.818 


5P ^ 5S 


0.638 


0.317 


2.605 



Table 8: The ratios among the radiative decay widths of the P-wave to the 5'-wave bottomonium 

states, where F = y(x2^~^t:i) ' -'^(^2 — >■ ^l) = rl^^Z^x-yj ^ ^^'^ ^^"^ values in the bracket come 
from the screened potential model |15j 



F 


X2 T27 


X2 T7 


Xi Tf27 


Xi T7 


Xo T7 


hb 7727 


2P ID 


0.156 [0.185] 


0.016 [0.013] 


0.527 [0.722] 


0.281 [0.268] 


0.561 [0.591] 


1.043 [2.371] 


3P ^ ID 


0.172 


0.013 


0.781 


0.295 


0.992 


1.163 


3P ^ 2D 


0.202 [0.190] 


0.016 [0.015] 


0.685 [0.635] 


0.279 [0.241] 


0.494 [0.471] 


1.046 [2.388] 


4P ^ ID 


0.175 


0.013 


0.815 


0.293 


1.059 


1.156 


4P ^ 2D 


0.185 


0.013 


0.815 


0.287 


0.941 


1.117 


4P ^ 3D 


0.204 


0.016 


0.620 


0.259 


0.433 


0.927 


5P ^ ID 


0.176 


0.012 


0.834 


0.294 


1.101 


1.164 


5P ^ 2D 


0.183 


0.013 


0.844 


0.291 


1.050 


1.144 


5P ^ 3D 


0.186 


0.013 


0.804 


0.285 


0.940 


1.106 


5P ^ 4D 


0.207 


0.016 


0.596 


0.245 


0.426 


0.903 



Table 9: The ratios among the radiative decay widths of the P-wave to the D-wave bottomonium 
states, where F = ^p7 — E^v^ — r, F(x2 — >■ T37) = = 1, and the values in the bracket come 

from the screened potential model [15] 
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i 


^2 — > X27 


2 Xi7 


T ^ X27 


Tf" ^ Xi7 


Tf" Xo7 


772 — > /l67 




n ofifi m 9/1 ni 

U.ZDD [U.Z4UJ 


n QQ/1 m snsi 
U.yy4 [U.oUoJ 


n 09/1 Fn 091^1 

U.UZ4 [U.UZOJ 


n /iKn rn /i9ni 

U.4DU [U.4ZUJ 


n 88"? rn k891 

U.ooo [U.DozJ 


1 19^ [1 /lonl 
i.izo [i.4yuj 




n 0/1 ^ rn 1 qoi 

U.Z40 [U.iozJ 


n Qi /I [1 1 of^i 
u.oi4 [i.iyoj 


n nof^ rn ni qi 

U.UzO [U.UioJ 


n /I QQ rn c;ni 1 

U.4oo [U.DUiJ 


n ^^oQ ri ^^^91 

U.oyo [i.DDzJ 


1 n^^n [1 0^71 
i.uou [i.yo<j 


op) V op 


U.Zoo [U.Z4UJ 


n S7/1 m 7R1 1 

U.o / 4 [U. < DiJ 


n noi rn 091^1 

U.UZO [U.UzOJ 


n /i/in rn "^ooi 
U.44U [u.oyyj 


n 8/1 1 rn 1^071 


1 1 n/1 [1 c;q71 


qp\ V 1 p 


n 0/1 7 


u. 1 yz 


n 097 

u.uz r 


n /1Q1 

U.4oi 


U.D40 


1 nqo 
i.uoy 


qp\ V OP 


n 0/1 Q 


n 700 
u. 1 yy 


U.UzD 


n /i9f^ 

U.4ZD 


U.D04 


1 n/1 9 

l.U4z 


3D ^ 3P 


0.230 


0.880 


0.022 


0.434 


0.877 


1.085 


4D ^ IP 


0.248 


0.781 


0.027 


0.429 


0.624 


1.029 


4D ^ 2P 


0.246 


0.781 


0.027 


0.425 


0.627 


1.027 


4D ^ 3P 


0.242 


0.801 


0.026 


0.427 


0.679 


1.033 


4D ^ 4P 


0.226 


0.953 


0.022 


0.471 


0.999 


1.185 


5D ^ IP 


0.249 


0.775 


0.027 


0.426 


0.610 


1.022 


5D ^ 2P 


0.247 


0.773 


0.027 


0.423 


0.608 


1.019 


5D ^ 3P 


0.245 


0.782 


0.026 


0.423 


0.633 


1.021 


5D ^ 4P 


0.240 


0.825 


0.025 


0.437 


0.707 


1.069 


5D ^ 5P 


0.222 


0.986 


0.021 


0.479 


1.006 


1.214 



Table 10: The ratios among the radiative decay widths of the I?- wave to the P-wave bottomonium 

states, where T = r, FfTa — >■ Y21) = tw^-7^^=1, and the values in the bracket come 

from the screened potential model [TS] 
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